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Abstract 



We present experimental results which demonstrate that nuclear magnetic 
resonance spectroscopy is capable of efficiently emulating many of the capabili- 
ties of quantum computers, including unitary evolution and coherent superpo- 
sitions, but without attendant wave-function collapse. This emulation is made 
possible by two facts. 

The first is that the spin active nuclei in each molecule of a liquid sample 
are largely isolated from the spins in all other molecules, so that each molecule 
is effectively an independent quantum computer. The second is the existence of 
a manifold of statistical spin states, called pseudo-pure states, whose transfor- 
mation properties are identical to those of true pure states. These facts enable 
us to operate on coherent superpositions over the spins in each molecule using 
full quantum parallelism, and to combine the results into deterministic macro- 
scopic observables via thermodynamic averaging. We call a device based on 
these principles an ensemble quantum computer. 

Our results show that it is indeed possible to prepare a pseudo-pure state 
in a macroscopic liquid sample under ambient conditions, to transform it into a 
coherent superposition, to apply elementary quantum logic gates to this super- 
position, and to convert it into the equivalent of an entangled state. Specifically, 
we have: 

• Implemented the quantum XOR gate in two different ways, one using 
Pound- Over hauser double resonance, and the other using a spin-coherence 
double resonance pulse sequence. 

• Demonstrated that the square root of the Pound-Overhauscr XOR corre- 
sponds to a conditional rotation, thus obtaining a universal set of gates. 

• Devised a spin-coherence implementation of the Toffoli gate, and con- 
firmed that it transforms the equilibrium state of a four-spin system as 
expected. 

• Used standard gradient-pulse techniques in NMR to equalize all but one of 

the populations in a two-spin system, so obtaining the pseudo-pure state 

that corresponds to |00). 

• Validated that one can identify which basic pseudo-pure state is present 
by transforming it into one-spin superpositions, whose associated spectra 
jointly characterize the state. 

• Applied the spin- coherence XOR gate to a one-spin superposition to create 
an entangled state, and confirmed its existence by detecting the associated 
double-quantum coherence via gradient-echo methods. 



2 



1 Introduction 



The theory of quantum computing is advancing at a rate that vastly outstrips 
its experimental realization (for accounts, see jl], ^, |^). Most attempts to imple- 
ment a quantum computer have utilized submicroscopic assemblies of quantum 
spins, which are difficult to prepare, isolate, manipulate and observe. A "ho- 
mologous" system that exhibits many of the same properties, but is easier to 
work with, would clearly be very useful both as a means of testing the the- 
oretical predictions, and exploring implementation issues like error correction. 
Such a system is provided by weakly polarized macroscopic ensembles of spins, 
which are readily manipulated and observed by nuclear magnetic resonance 
spectroscopy, or NMR. 

The spins of a molecule in solution are largely isolated from their surround- 
ings by simple surface-to- volume considerations, and from the spins in neighbor- 
ing molecules by diffusional motion, which averages their dipole-dipole coupling 
to a second-order effect . This fact enables us to work with a reduced density 
matrix ^ of size 2", where n is the number of spin ^ nuclei in the molecule, 
rather than 2^ where N is the total number of such spins in the sample . It 
is also customary in NMR spectroscopy to shift the reduced density matrix by 
subtraction of its mean trace, since only the traceless part undergoes unitary 
evolution, and to scale it to have integral elements Q. In the next paragraph, 
we define a manifold of statistical spin states with a reduced density matrix 
whose traceless part is proportional to the traceless part of the usual density 
matrix of a pure state. 

Henceforth, whenever we use the term "density matrix" , we mean "reduced, 
shifted and scaled density matrix" unless otherwise stated. When such a den- 
sity matrix has rank equal to one (after adding an appropriate multiple of the 
unit matrix to it), it can be factored into a dyadic product of the coordinates 
of a "spinor" and its conjugate versus the usual Iz basis, and this factorization 
is unique up to an overall phase factor. This mapping between spinor coordi- 
nates and density matrices that can be shifted to a signature of [±1,0, .. . ,0] 
is covariant, in the sense that if we apply a unitary matrix to the spinor's 
coordinates, the corresponding density matrix transforms by conjugation with 
the same unitary matrix. As a result, we can regard such a density matrix 
as a kind of spinor, and perform essentially arbitrary unitary transformations 
on it via NMR spectroscopy, thereby "emulating" a quantum computer. We 
shall call the states described by density matrices with 2" — 1 equal eigenvalues 
"pseudo-pure" states, and the corresponding spinors "pseudo-spinors" . 

Of course, some things are lost in translation. For example, the density ma- 
trix is not changed on rotation by 27r, although spinors change sign. Since these 
sign changes cannot easily be observed, this seems to be of little consequence for 
quantum computing. More important is the fact that the "coherence" observed 
by NMR spectroscopy is always an ensemble average over an astronomical num- 
ber of microscopic quantum systems. As a consequence, the NMR spectrum of 
a pseudo-pure state yields the expectation values of certain observables relative 
to the corresponding pseudo-spinor, rather than a random eigenvalue of one 
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of them. In particular, wave function collapse does not occur. A variety of 
other more easily controlled "filtering" mechanisms are available in NMR spec- 
troscopy, however, and we have shown that for most computational purposes the 
ability to measure expectation values directly is actually a great advantage Q . 
NMR experiments on liquid samples possess a number of other highly desirable 
features as well; in particular, the decoherence times are typically on the order 
of seconds. 

NMR spectroscopy in fact provides a means of building a nonconventional 
computer that can be programmed much like a quantum computer, but is much 
easier to implement on at least a limited scale. In some respects, this approach 
also resembles DNA computing, in that it can use the parallelism inherent in 
ensembles of molecules to efficiently count the number of solutions to combina- 
torial problems, trading an exponential growth in the time required against an 
exponential growth in the sample size. More generally, we have called a compu- 
tational device that operates by running a large number of quantum computers 
on coherent superpositions, and then estimates the expectation values of observ- 
ables by summing them over all the quantum computers, an ensemble quantum 
computer. A detailed introduction to the theory of such machines may be found 
in 1^; this paper will describe how basic quantum logic gates can be imple- 
mented via NMR spectroscopy, and present experimental results to validate our 
claims. After the majority of these results had been obtained we learned of 
a similar approach proposed by other researchers ^ . 

2 Basic techniques from NMR 

This section introduces the basic techniques from NMR spectroscopy that are 
needed for this paper, and in the process defines the notation it uses (for more 
complete introductions, see e.g. [p|, p|, p^). 

Let us consider the simplest nontrivial case, which contains all the essential 
ingredients of solution NMR spectroscopy. This is a liquid consisting of identical 
molecules each containing exactly two coupled spin ^ nuclei of the same isotope 
(throughout this paper, ^H). The dipolar coupling between the spins is averaged 
to zero by the rotational motion of the molecules in the liquid, and hence the 
coupling in this case is the so-called scalar coupling, which is mediated by 
electron correlation in the chemical bonds linking the atoms. It will simplify 
our presentation if we assume weak coupling, i.e. that the coupling constant 
Ji2 is small compared to the difference jcji — (JJ2I m the resonance frequencies 
of the two spins. With the convention that the magnetic field is along the z- 
axis, the Hamiltonian of this system is H = wilj -I- uj2'i-t + 27r J12 I^I^, where 
(fc — 1, 2) are the usual matrices for the z-component of the angular momentum 
of the spins. Because the energy level differences are small compared to kT at 
room temperature, the equilibrium density matrix of this system is given to an 
excellent approximation by Exp(—H/fcT) w 1 — H/fcT. On taking account of 
the fact that | J12I <ti ui ~ uj2 , removing the trace and scaling, the equilibrium 
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density matrix becomes 
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(where the longitudinal spin states which label the rows and columns of this 
matrix are shown along its right-hand side). 

Rather than writing them out explicitly, NMR spectroscopists typically rep- 
resent their density matrices as linear combinations of "product operators", i.e. 
products of the usual angular momentum operators ij, 1^, (as in Eq. (^) 
. This makes it very easy to describe the unitary transformations effected 
by applying RF (radio-frequency) pulses to the sample. For example, a "soft" 
pulse whose frequency range spans the resonance frequency of only the first spin, 
and which imparts an energy sufficient to rotate that spin by an angle (/? about 
the y-axis (in a frame rotating with the carrier of the receiver Q ) is 

lii^cos(^)li+sin((p)li . (2) 



Thus when a "soft" (i.e. spin-selective) [tt/2] y-pulse is applied to the first spin of 
a two-spin system at equilibrium, we obtain l!^-f I^, while a "hard" (nonselective) 
[it/2] y-pulse yields l^+l^.. 

The density matrix evolves according to the time-dependent unitary trans- 
formation 1^ 

^ ^ Exp( i t H) * Exp(i t H) . (3) 

Since all three terms of the Hamiltonian commute, the above propagator fac- 
tors into a product of the chemical shift and scalar coupling propagators. The 
chemical shift propagator for the first spin can be expanded as 

Exp(iiwil^) ^ cos(cJii/2)l + 2isin(wit/2)I^ . (4) 

An exercise in the Pauli matrix algebra then shows that 

Exp(-ittJiI^) I;^ Exp(iiwil^) = cos(wit)li -f sin(u;ii) . (5) 

Altogether, we obtain: 

> cos(a;it) l]. -\- sin(a;ii) 
untV^ cos(t^it) ll - sin(wii) ll (6) 



This propagator has no effect on terms involving only the second spin, which 
evolve analogously under their own chemical shift Hamiltonian. 
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The scalar coupling propagator, on the other hand, is 

Exp(27ri Jiatl^I^) = cos(7r Ji2V2)l + 4i sin(7r Ji2i/2) I^I^ . (7) 
A similar calculation shows that it transforms the one-spin operators as follows: 



I 27rJi2tI^lj 



1 27TJr2tllll 



cos(7rJi2t) ll + 2 sin(7r Ji2i) I^I^ 



cos(7rJi2t) I,y - 2 sin(7rJi2t) I^I^ 



(8) 



with analogous expressions for the terms I^, and I^. 

Physically, the above expressions describe the precession of the transverse 
magnetization about the applied field, which generates a rotating magnetic mo- 
ment in the xy-plane. The complex-valued signal induced in the receiver is 
calculated by taking the trace of the product of this time-dependent density 
matrix with 
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where = + il^ as usual (see e.g. [|l3|). Since this matrix contains only 
four nonzero elements, the spectrum contains direct information on only four of 
the ten independent elements of the density matrix, namely 4'i2 = ^'2i,^'i3 = 
^311^24 = ^42 sind \l/34 = '^X^,- These elements are called single- quantum 
coherences^ because they connect pairs of states related by single spin flips (cf. 
Fig.|). 

For example, the signal due to the chemical shift precession of lij, alone is 



tr(l^ (cos(wit)I;J, -|- sin((jjit)lj)) — cos(ti;it) -|- i sin(ajit) . 
These two terms are modulated by the scalar coupling evolution to 

cos(wit) tv{l\ {cos{-KJi2t)ll + 2 sin(7r Ji2i) ijl^)) , 

= cos{uJit) cos(7r Ji2t) 
sin(a;i t) tv{l\ {cos{-KJi2t)ll - 2 sin(7rJi20 lllD) 

= sin(ajit) cos(7rJi2t) . 

The total signal due to the first spin is thus 

cos{u!it) cos{TTJi2t) + sin(ijjii) cos(7rJi2t) 
= i(exp(i(wi - 7rJi2)t) + exp(i(tJi + 7rJi2)i)) , 



(10) 



(11) 



(12) 



which shows that the real part of the Fourier transform consists of a pair of 
peaks centered on the resonance frequency of the spin and separated by the 
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Figure 1: The experimental NMR spectrum of liquid 2, 3-dibromo-tliiophene 
(see below) in a 9.4 Tesla field, at which the two doublets are separated by 130 
Hz. This spectrum was collected by applying a nonselective [tt/2] y-pulse to 
the equilibrium state and Fourier transforming the resulting signal (see text). 
Note that frequency (the horizontal axis) increases from right-to-left in NMR 
spectra. 



coupling constant; this is called an in-phase doublet, meaning that both peaks 
have the same sign. Thus if the magnetization due to both spins is rotated into 
the transverse plane by a hard [tt/2] y-pulse, one obtains a spectrum containing 
a pair of doublets. This is shown in Fig. ^ using the molecule 2, 3-dibromo- 
thiophene shown in Fig. ^. 

Once the magnetization due to a spin has been placed in the transverse 
plane, a soft [tt] pulse in the middle of an evolution period t may be used to 
refocus its chemical shift evolution during that period, as follows: 

II cos{u;it/2)ll + sin(c^ii/2)li 



LUitll/2 



(13) 



- cos{LJit/2)ll + sin(cjiV2)li 

— COs((jJi t/2)(cos(cJit/2)li + sm{LUit/2)ll) 

+ sm{uJit/2){cos{uJit/2)ll - sin(wit/2)I^) 

Standard trigonometric identities show that this last expression is simply — 
Thus, up to an inconsequential overall phase factor, we have cancelled the effect 
of chemical shift evolution. 

A very similar calculation shows that such a soft [tt] y-pulse also cancels the 
scalar coupling evolution, so that only the chemical shift evolution of the other 
spin occurs during the interval, i.e. 

X ' X 

" II + cos(^2t)I^ + sm{uj2t)ll . 

We shall denote such a rotation of the second spin about the z-axis by [oJ2t]1- 
Note, however, that [tt]1 is more easily obtained as [tt]^ — [tt]^ (a [n]y pulse 
immediately followed by a [tt]^). 
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Figure 2: The chemical structure of 2,3-dibromo-thiophcnc. This molecule is 
a liquid at room temperature, and contains two inequivalent, spin i hydrogen 
atoms with a coupling constant of 6.0 Hz. 

A hard [tt] y-pulse applied in the middle of a period, on the other hand, 
cancels the chemical shift evolution of both spins while allowing their scalar 
coupling to evolve. This is because the pulse just changes the sign of the density 
matrix half-way through the period, i.e. 

ji ^J^^t-^l^l^ cos{nJi2t/2)ll + 2 sin(7rJi2i/2)lJ,l' 

- cos(7r Ji2t/2)li - 2 sin(^Ji2t/2) I^I^ (15) 



' ' — cos(7r Ji2t)I^ — 2 sin(7r Ji2t) ijl^ , 

thereby yielding the negative of what one would have obtained with no [tt] pulse 
(as shown). We will denote such a scalar coupling evolution, with no chemical 
shift evolution, by [t\. 

A gradient pulse produces a transient field inhomogeneity along the 2;-axis, 
which has the Hamiltonian 

Hgrad = jr-VB, (16) 

where r is the position vector within the sample. Although the microscopic evo- 
lution remains unitary (of course!), a gradient pulse dephases the macroscopic 
transverse magnetization due to the off-diagonal elements of the density matrix. 
The decay rate of each element is proportional to the difference in the number of 
"up" spins between the corresponding pair of states, which is commonly called 
the coherence order. The net effect is to zero the off-diagonal elements of the 
density matrix whose coherence order is nonzero, thus effecting a projection of 
the state. 

This makes it possible to destroy (or more precisely, render unobservable) 
the magnetization due to selected spins with a combination of soft and gradient 
pulses, e.g. 

Il + Il ^ Il+Il II. (17) 
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tt 








Figure 3: The energy level diagram for a coupled two-spin system, where 
the resonance frequencies of the two spins are vi and i>2, while the coupling 
constant is J12. The single-quantum transitions allowed by the selection rules 
for angular momentum are indicated by two-headed arrows, together with the 
energy associated with each. 



A further [tt / 2] y-pulse transforms the remaining term to I]. , which is converted 
by a scalar coupling evolution of [1/(2 J12)] to 21jl^, and thereafter evolves as 

cos(cJii) (2 cos(7r Ji2t) - sin(7r Ji2t) I;^) 

19 1 \ \^°J 

- sm{ujit) (2 cos(7rJi2t) I^If + sin(7r Ji2i) 1^) 
Taking the trace product of this with + I'^ yields the signal 
— cosiujit) sin(7rJi2i) — sinfo^it) sin(7r Ji2i) 

(19) 

= 2 ('5^P(i('^i ^ ^Ji2)t) - exp(i(tJi + 7rJi2)i)) . 

Thus if one collects a spectrum starting after the [1/(2 J12)] evolution, one ob- 
tains an anti-phase doublet consisting of two peaks of opposite sign. 

Alternatively, one can apply a selective [— 7r/2]^ pulse to the second spin of 
2IyI^, obtaining the correlated state 21jl^, or a [n/2]y pulse to obtain 21^1^. 
These contain no single- quantum coherences, and hence produce no observable 
magnetization. 



3 The Pound- Over hauser XOR and conditional 
rotation 

The most obvious way to implement the quantum XOR (or controlled-NOT) 
gate is to use a pulse that is selective for just one component of a doublet; this 
effects a population transfer similar to the original ENDOR experiment, and 
constitutes an example of Pound- Overhauser double resonance Specifically, 
if we apply a [tt] pulse about the y-axis to the transitions 1 <-^- 3 and 2 <-!■ 3 
(see Fig. H), we obtain an XOR gate with the output on the first and second 
spins, respectively. These pulses will be denoted by [7r]+'^ (fc = 1,2). If |ei,e2) 
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Table 1: Simplified "stick" spectra for each of the three diagonal product op- 
erators that would be observed after selective [7r/2]^ observation pulses on the 
/c-th spin, before (initial) and after (final) applying an XOR with its output on 
the second spin (see text). 



denotes the spinor of a single state (ei,e2 G {0,1}), these pulses effect the 
unitary transformations of the density matrix of the corresponding pure state 
given by 

Up''o|ei,e2)(ei,e2|Upo = |ei ® £2, e2>(ei ® £2, £2! 

(20) 

Up''o|ei,e2)(ei,e2|U|,o = ki, ei © £2) (ei, ei © £2! , 
where "©" denotes the boolean XOR operation, and 
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(21) 



By applying [tt] pulses to the transitions <-> 1 and <-> 2, we obtain the boolean 
operations \ei © £2, £2) and |£i, £1 © £2) (where the overbar denotes the NOT of 
the corresponding qubit); these pulses will be denoted by [tt]"*"' (k = 1,2). 

The necessary selectivity can be obtained with a long sinc-modulated pulse, 
whose Fourier transform has a square-wave envelope occupying just the width 
of a single peak. Based on the matrices in Eq. (|2l] ) above, the result of applying 
such a pulse to the first spin (i.e. the left-most peak in its doublet) of a two-spin 
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Figure 4: The spectra obtained from [7r/2]^ readout pulses on the first (above) 
and second (below) spins, following a [7r]+^ pulse applied to the equilibrium 
state (see text). 

system at equilibrium should be 

|00)(00|-|ll)(ll|=li+l2 

r ,+1 (22) 

^ mm-\01){Ql\^2llll+ll. 

This expectation can be confirmed by collecting spectra following ordinary soft 
readout pulses on each spin. The spectra expected from the individual product 
operators are shown in a diagrammatic "stick" form in Table 0, and the spectrum 
which results from any sum of these operators will be the point-by-point sum of 
the spectra from the individual terms in the sum. For example, a readout pulse 
on the second spin yields 

2lil^+I^ ^ 2l\ll + ll (23) 

Due to interference between the in-phase and anti-phase signals, the resulting 
spectrum contains only a single peak, with twice the intensity of the peaks one 
gets from a readout on the first spin, as shown in Fig. |[ 

Similarly, the result of applying a Pound-Overhauser XOR to the second 
spin is + 21^1^. The effects of a Pound-Overhauser XOR on all the two-spin 
product operators are shown in Table ^. Note that the matrices in eq. (|l]) differ 
in the sign of one element from the matrices for the XOR usually encountered 
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Table 2: The effect of the [ttJ^ ^ pulse on ah the product operators of a two- 
spin system. The two factors of the product operator on which the XOR acts 
are shown in the left column and top row, while the result is shown in the 
corresponding table cell. Thus for example, a [7r]+^ pulse converts the product 
operator to 2 1;^,!^ . 



in quantum computing, e.g. 
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(24) 



This sign difference has no effect on the results of applying the gate to the 
diagonal product operators (i.e. single states), but the results of applying the 
Pound-Overhauser XOR to a superposition may differ by phase factors from the 
results obtained with the quantum computing XOR. 

To see how to get the same phase factors, we compute the infinitesimal 
generators of the two XOR's, i.e. 



Upo = Exp(i7rOpo) where Gpo 



and 





/o 








'\ 


1 











— i 


2 
















vo 


i 



















1 





1 





-1 


2 
















^0 


-1 





1/ 



(25) 



Uqc = Exp(i7re^c) where Gjjc = 



These generators may be expressed in terms of product operators as 
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Table 3: The effect of the [XOR]qq pulse sequence on all the product operators 
of a two-spin system (cf. Table ||). 



in both of which all the terms commute with one another. The term jl just leads 
to an overall phase shift, while the term — 5I2 can be cancelled with a [7r/2]2 
rotation. It follows that the unitary transformation Uqq can be obtained with 



1+1 



pulse together with a [— 7r/2]2 rotation (in either order). The effect 
of this XOR]qq pulse sequence on the two-spin product operators is shown in 
Table |. Comparison with Table ^ shows that the two types of XOR operations 
differ by swaps of ^1 with in those product operators that contain either of 



these factors. 



^1 . 

and Ir with I- 



in all the rest. 



The XOR gate together with arbitrary one-bit rotations, which one can im- 
plement via soft pulses and free precession, constitute a universal set of quantum 
logic gates. It is nevertheless of interest to observe that the Pound-Ovcrhauser 
XOR gate generalizes directly to a conditional rotation. For example, the square 
root of the XOR is obtained by applying a [7r/2]+^ pulse (i.e. the same pulse 
needed to get the XOR, but for only half as long). For the output on the first 
bit, this leads to the matrix 
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(29) 



To get a unitary transformation whose matrix is that of the square root of the 
quantum computing XOR, i.e. 



1 + i 
V2 



/I -i 
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(30) 



one need only apply a [— 7r/2]+^ pulse together with a [— 7r/4]2 pulse (in either 
order) . 



Based on the matrix in Eq. (29), the effect of a [7r/2]+^ pulse on the equi- 
librium state should be: 
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(31) 
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Figure 5: Validation of a Pound-Overhauser conditional rotation. The upper 
spectrum was collected immediately after applying the [7r/2]+"'^ pulse to the equi- 
librium state, while the lower spectrum was collected following the application 
of an additional hard [7r/2]^^ pulse. 

This gives rise to a spectrum containing a single peak at the very same frequency 
that the pulse was tuned to. An additional hard readout pulse produces the 
state 

-^Il+llll + llll+'^ll + ll . (32) 

The corresponding spectra are shown in Fig. ^. The general formula for a 
conditional rotation in terms of product operators is: 

ll cos^(^/2)li+cos(v./2)sin(^/2)(li-2lil^) (33) 

+ 2 sin2(^/2)lil2 

4 The spin-coherence XOR and ToffoU gates 

Due to the degree of selectivity required and the necessity for direct coupling 
between the spins involved, the Pound-Overhauser XOR gate can be difficult 
to apply. We have therefore developed the following pulse sequence, which 
constitutes an example of spin- coherence double resonance [p!^, and promises 
to be more generally useful: 

[XORfsc ^ i^/ty - [1/(2 Ji2)] - k/2]J: (fc = 1, 2) (34) 
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Table 4: The effect of the [XOR]\q pulse sequence on all the product operators 
of a two-spin system (cf. Tables 2 & 3). 



To show that this does indeed effect the boolean XOR operation, we demonstrate 
that it has the same effect on the diagonal product operators as the Pound- 
Overhauser XOR above, e.g. for fc = 1: 



y2 1^/21 « j2 [l/(2Ji2)] 2 .2 

z z z z 

21111 ^21111 [V(2.M1 ^ JVgl^ I, 



(35) 



Table ^ shows the effect of this gate on all the two-spin product operators. 
Alternatively, we can just multiply together the matrices of the three individual 
steps, to obtain 
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(36) 



We shall call this the spin- coherence XOR gate. 

The validity of this XOR implementation is corroborated by the spectra 
shown in Fig. ^, which were obtained by applying this sequence to the equi- 
librium state followed by soft readout pulses. The infinitesimal generator of 
the product of the matrix in Eq. ( ^6| ) with the matrix of the usual quantum 
computing XOR from Eq. ( |2^ ) is 



Exp(i7rA) where A = il 



2^2 



il2 

2 2 



(37) 



Together with the fact that Uqq is self- inverse and commutes with Exp(i tt l2/2), 
this generator shows that the phases of the nonzero components can again be 
equalized by composing this pulse sequence with suitable z-rotations, namely 
[',/2]l-[XOR]lc-[-^/2]\. 

We have also developed a pulse sequence, analogous to the above [XOR]'^^ 
sequence, which transforms the longitudinal spin states of a three-spin system 
according to the truth table of the well-known Toffoli gate flS] . We call this the 
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Figure 6: The spectra obtained from readout pulses on the two spins, following 
a [XOR]gQ pulse sequence applied to the equilibrium state (see text). 



Toffoli pulse sequence: 

[TOF]lc^ [7r/2]^-[l/(4J)]-[7r/2]^-[l/(4J)] 
-[-7r/2]^-[l/(4J)]-[-7r/2]^ 

This pulse sequence assumes that the coupling constants Jke and Jkm have 
the same value J, which can always be arranged by inverting the spin whose 
coupling constant with k is greater using a soft [tt] pulse part way through each 
delay. This changes the sign of the effective coupling constant, so that the time- 
average coupling constant can be given any desired value between the original 
coupling constant and its negative. 

If we place the output on the first spin {k = 1), the matrix of the above 
Toffoli sequence can be shown to be: 



/1+i 























\ 





i 




























i 








































1 - 


i 














1-i 




























— i 




























— i 







I 








-1-i 














J 



(39) 



This is easily seen to act upon the density matrices associated with single states 
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Figure 7: The chemical structure of l-chloro-2-nitro-benzenc, which was used 
to vahdate the spin-coherence implementation of the Toffoli gate by applying 
the [TOFjgQ sequence to hydrogen atoms 1 — 3. The relevant coupling constants 
among these atoms are J12 = 8.0 Hz., J13 = 1.5 Hz., and J14 < 0.1 Hz., and 
J23 = 7.0 Hz. 



as 



Vgj.|ei,e2,e3)(ei,e2,e3|Vsc 
= |ei ® (£2 A €3), £2, e3)(ei ® {^2 A €3), €2, £3! 



(40) 



(where the "A" denotes the boolean AND operation). Unhke the XOR sequence, 
however, the Toffoli sequence does not simply permute the diagonal product 
operators; instead, it results in a sum of such diagonal operators, e.g. 



li 



k/2]i 



[1/(4J)1 



[1/(4J)] 



[1/(47)1 



[-T/2]i 



2'-x ^ '-y'-z ^ '-y'-z ^ '-x^z'-z 



2^1^ ^y^z + ^y^z + ^ I^I^I^ 



riT2T3 



ill _ ill 

2 z 2 X 



2jlj2j3 

X z z 



2 Jll2j3 

z z z 



2^ 2^x '^^x^z^z '^'^y^z^t 



2-^y 



Ill2 

y z 



2 z ' z z ' z z 



rlrS 



2jlj2j3 

z z z 



(41) 
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Similar calculations lead to the following complete list: 





[TOF 


Isc 
> 


2 z ^ ^z^z ^ ^z^z 


- ^iim 


ll 


[TOF 


Isc 
> 








[TOF 


Isc 
> 


^z 






[TOF 


Isc 
> 


Ijl 1 t1t2 t1t3 

2 z ^z^z ^z^z 


+ 21^1^13 




[TOF 


Isc 


t1t2-i-T1t3 

2 z ^z^z ^ ^z^z 


+ 2lil^l3 


21^1^ 


[TOF 


l^c 


Z Z 




4111^1^ 


[TOF 


Isc 


~ 'qJ'Z ^z^z ~l~ ^z 


+ 21^1^1; 



(42) 



Once again, these claims can be corroborated by using the matrix in Eq. ( p9D 
to predict the result of applying the [TOF]^^ sequence to the equilibrium state, 
and then verifying that the results are consistent with the spectra collected after 
appropriate readout pulses. These spectra are rather complicated, however, 
and hence we shall simply show the result of applying the readout pulse to the 
first spin. The molecule used for these experiments, l-chloro-2-nitro-benzene, 
is actually a four-spin system, thus enabling us to also demonstrate that we 
can apply quantum logic gates to subsets of spins, given sufficient frequency 
resolution. This spectrum is shown in Fig. ||, along with the spectra obtained 
by applying the same readout pulse to the four diagonal product operators 
(prepared by the gradient-pulse techniques described earlier) for comparison. 
The fact that just one of the four peaks due to the first spin has been inverted 
by the Toffoli sequence is a direct reflection of the fact that the population 
difference of the first spin has been inverted in just those molecules wherein the 
other two spins were "up" . 

In order to equalize the phase factors, we compute the infinitesimal generator 
of the product of the matrix Vg(^ for the [TOF]^^ gate shown in Eq. ( |39| ) with 
the desired matrix (consisting of ones at all the same nonzero locations): 

Exp(i7rA) = VjjcVgc where 

A = il-li-il2_ il3_ iiii2_ ijiis , Il2j3 ^ 

8 -^2: 4^z 4-^2: 2 z'-z 2 z'-z ' 2 z^z 

These are all z-rotations of one kind or another, and hence can all be readily 
obtained from the three-spin analogues of the implementations introduced above 
for a two-spin system (save possibly for the rotation whose generator is I^I^, 
which will require a relay sequence if these spins are not directly coupled p^). 

The Toffoli gate is known to be universal for classical computation, and 
hence the existence of this pulse sequence shows that any boolean function can 
be implemented in NMR via soft pulses selective for single spins. A Pound- 
Overhauser implementation of the Toffoli gate is also possible, but we have not 
yet done the experiments to demonstrate this. 
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Figure 8: The spectra of l-chloro-2-nitro-benzene obtained following a [7r/2]^ 
readout pulse to the four diagonal product operators I^, 21^1^, 21^1^ and 
21^1^1^, as well as to the result of applying the [TOFjgg pulse sequence to 
the equilibrium state (see text). 
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5 Superpositions of pseudo-spinors 



If one generates all three of the states I^, and 2I]Jl simuhaneously in the 
same sample, one obtains the density matrix 



/ ^ 

' 2 



This matrix shifts to 



21^1? 








\ 






-y 



(44) 



/ 2 \ 





\ y 



(45) 



which can in turn be factored into the product of a pseudo-spinor with its 
conjugate, 



2|00)(00| 



/1\ 







(1 ) 



(46) 



Therefore, the expression + + 21^1^ in Eq. (^) represents a pseudo-pure 
state in product operator notation. 

The following RF and gradient pulse sequence transforms the equilibrium 
state of a two-spin system into this pseudo-pure state: 



[grad]j 



11+11/2 ~llV3/2 



[l/(2./l2)] 



[grad]j 



li/V2 + l2/2-li/x/2 
1I/V2 + II/2 + V2II1I 
ll/2 + ll/2-ll/2 + llll 
\\I2 + \II2 + \\\1 



(47) 



jll2 

z z 



Subsequently, any of the four basic pseudo-pure states (with a density matrix 
shifting to a matrix with a single nonzero element on the diagonal) can be 
obtained by means of soft [tt] pulses (which perform the NOT operation on the 
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spins). In terms of product operators, these can be written as: 



2|00)(00|-il = Il+ 1^ + 2 III, 



2|01)(01|-il = Il-Il^2llll 
2|10)(10|-il = -II+II-21I1I 
2|ll)(ll|-il = -II-II + 2II1I 

The factor of 1/2 difference between the last hnc of Eq. ( ^7|) and the first 
fine in Eq. (^ ) is due to a 25% loss of total polarization during the gradient 
pulses. Note that by changing the sign of cither of the two [tt/A] rotations in 
Eq. (^tI), we can generate the state I?; + — 2 This is the negative of last 
state in Eq. (^8|), which means the polarization of the sample (i.e. the sign of the 
nonzero element after shifting the density matrix) is inverted between the two 
states. The negatives of all four states in Eq. (|4|) can likewise be interconverted 
by soft [tt] pulses, but none of these states can be converted to its negative by 
means of RF pulses. Whether we use the four basis states in Eq. ( ^ ) or their 
negatives is irrelevant for computational purposes. 

To validate that with the use of a [7r/4]J instead of a [— 7r/4]^ pulse in 
Eq. (^7|), we have actually created the expected state 12+12 — 21^1^, we consider 
the spectra that are obtained after applying soft readout pulses to each of the 
two spins, i.e. 

11+1^-2111^ ^ ll+ll-2llll (49) 



and 



ll+ll-2llll ^ ll+ll-2l\ll. (50) 



This creates the difference of an in-phase with an anti-phase state in both cases, 
so that the associated spectra consist of single peaks at the leftmost position of 
the corresponding doublet. These spectra are shown in Fig. |^. 

The density matrices generated by these readout pulses can be written in 
terms of the associated pseudo-spinors as 

il-(|0)-|l))|l)((0|-(l|)(l| ^ ll + ll-2llll 
il-|l)(|0)-|l))(l|((0|-(l|) = Ii+l^-2lll2 . 



In other words, each readout pulse puts the corresponding spin into a super- 
position over its longitudinal states. We can create a superposition over both 
spins by performing the two readout pulses in rapid succession, or with a single 
"hard" pulse (which takes only a small fraction of the time required for a soft 
pulse), i.e. 

il_2|ll)(ll|i^li +1^^2111^ (52) 
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Figure 9: Two experimental NMR spectra of 2, 3-dibromo-thiophene, which 
corroborate the creation of the pseudo-pure state + — 2l],I^ via readout 
pulses selective for the first (above) and second (below) spins (see text). 



Since 2 does not contribute to the signal, the spectrum in this case consists 
of a pair of in-phase doublets, exactly like that shown in Fig. |l]. 

In NMR computing, superpositions over the basic pseudo-pure states are 
always associated with nonzero off-diagonal coherences in the density matrix. 
The product operator above, for example, represents the matrix 



Ii + 12 - 



21^1^ 

X X 



= i(l-(|0)-|l))(|0)-|l)) 
((0|-(1|)((0|-(1|)) 
■ 1 1 - 1 \ 







1 



1 



V 



1 



1 



(53) 



Since the corresponding pseudo-spinor can be factored into a product of one-spin 
pseudo-spinors (as shown), it represents an "unentangled" state. 

The fact that we have successfully prepared a pseudo-pure state, and our 
spin-coherence XOR gate, can both be further corroborated by putting them 
together. The resulting spectra are shown in Fig. 
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We have also confirmed 
that the superposition created by applying a soft [n/2]y pulse to the ll + + 
21^1^ state, namely (|0) -I- |1))|0), is changed by only inconsequential phase 
factors on applying the [XOR]^q sequence to it. These phase changes can 
nevertheless be made visible by applying a further [7r/2]^ readout pulse to the 
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Figure 10: The eflFect of applying the [XOR]^^ pulse sequence to the pseudo- 
pure state II +Il — 2IIII to get —I] + + 21^1^. This is clearly seen from 
the change in sign of the first in-phase doublet following a hard [n/2]y^ pulse, 
before (above) and after (below) the XOR. 

result, as shown in Fig. ^ and described in its caption. Finally, we consider 
the result of applying the [XOi?]g(^ gate, with its output on the other spin, to 
this same one-spin superposition. As is well-known, this creates an "entangled" 
state 



II 



[XOB^fsc 



(54) 



or equivalently: 



+ (|00) + |10)) [xoR]l, 
m\ + {10\)) 



-1 + (|00) + |11)) 

((oon-(iil)) 



(55) 



In this case the transverse magnetization corresponds to a double- quantum 
coherence. 



2 '^x^x 



2T^r — 



/ 





VI 



1\ 





/ 



(56) 
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Figure 11: The result of applying the [XOR]gQ sequence to the one-spin 
superposition created by applying a [7r/2]J pulse to the pseudo-pure state -f 
II + 2IIII, to get + 21jl^ (above). This is then further confirmed with 
an additional [7r/2]^ readout pulse to get ly + + 21^1^, which consists of 
two in-phase doublets that are 90° out-of-phase with each other, so that one 
doublet appears as a pair of dispersive peaks whenever the other is phased to 
be an absorptive pair (below). 

which precesses at twice the rate of the single-quantum terms, but produces no 
observable magnetization. The most direct way to confirm the creation of this 
double-quantum coherence is via a gradient echo experiment. In this technique, 
one apphes a -l-z-gradient for a period of l/(2Ji2), followed by a soft [— 7r/2]^ 
pulse, and then an inverse — z-gradient. The double-quantum coherence de- 
phases at twice the rate of a single-quantum coherence, but even though no 
macroscopic magnetization remains at the end of the -|-z-gradient, the micro- 
scopic coherence 2I;J,I^ is still converted to 2Ii^I^ by the 7r/2]^ pulse. This 
term then rephases under the — z-gradient at the usual rate for a single-quantum 
coherence, while at the same time evolving under coupling to 1^. This in turn 
results in observable transverse magnetization, which reaches its maximum after 
of period of 1/Ji2- The double-quantum echo appears after twice the time of 
the echo due to the rephasing of the (residual!) single-quantum coherence. If 
we collect the data for a spectrum starting right after the [— 7r/2]^ pulse, we 
observe the expected anti-phase doublet (Fig. [l^ ). 

For completeness. Table || gives the coefficients of the seven diagonal product 
operators in the eight basic pseudo-pure states of a three-spin system. A pulse 
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Figure 12: Time-domain echo of the double-quantum coherence character- 
istic of an "entangled" pseudo-spinor (left), together with the corresponding 
frequency-domain spectrum (right). The first echo, due to residual single- 
quantum coherence, occurs in half the time required for the double- quantum 
echo to form (see text). 



sequence and gradient sequence that generates the |000) state is: 

jl -L l2 J. i3 

[57r/12]i-[7r/3]J-[grad], 



[V4]'-[l/(2Ji2)]-[-7r/4]^-[grad]. 



[7r/4]^-[l/(2J)]-[7r/4]?-[grad], 



k/4]J-[l/(4J)]-[^/4]3-[gradl. 



4 z ^ 2 z ~ -^z 



4 z ^ 4 z ' 2 z' 2 z'-z ' 



(57) 



i'-z ^ i'-z ^ i'-z ^ 2 z'-z ^ 
2 z'-z ~ 2 z'-z ~ '-z'-z'-z 

In this sequence, [l/(2Ji2)] stands for a coupling evolution period with the 
coiipling constants J13 and J23 averaged to zero by a [tt] pulse on the third spin, 
while [1/(2 J)] and [1/(4 J)] are evolution periods with J12 averaged to zero and 
■^13, •^23 averaged to the same value J. 



6 Conclusions 

We have demonstrated that nuclear magnetic resonance spectroscopy provides 
an experimentally accessible paradigm for quantum computing. In particular, 
the results given in this paper include the first physical implementations of 
all the basic quantum logic gates, including the XOR and Toffoli gates, which 
have up to now been largely theoretical constructions. We have further shown 
that one can actually prepare a macroscopic ensemble of weakly polarized spin 
systems in a pseudo-pure state, which can be described by a spinor just like a 
true pure state, and confirmed that one can put these states into the equivalent 
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Spinor 


ll 






21^1? 


2I'I^ 


2I'I' 




|000) 


1 


1 


1 


1 


1 


1 


1 


|001) 


1 


1 


-1 


1 


-1 


-1 


-1 


|010) 


1 


-1 


1 


-1 


1 


-1 


-1 


|011) 


1 


-1 


-1 


-1 


-1 


1 


1 


|100) 


-1 


1 


1 


-1 


-1 


1 


-1 


|101) 


-1 


1 


-1 


-1 


1 


-1 


1 


|110) 


-1 


-1 


1 


1 


-1 


-1 


1 


1111) 


-1 


-1 


-1 


1 


1 


1 


-1 



Table 5: Table of coefficients of the seven diagonal product operators corre- 
sponding to each basic pseudo-spinor for a three-spin system. 



of entangled superpositions. Finally, we have shown how the state of the system 
can be efficiently determined from spectra collected following suitable readout 
pulses. 

These operations provide essentially all the ingredients needed to efficiently 
emulate a true quantum computer by NMR spectroscopy. Nevertheless, as has 
recently been forcefully pointed out []l6| , the population difference for any single 
spin that results from preparing a pseudo-pure state by averaging over unitary 
transformations of the thermal equilibrium state falls off as at least n/(2" — 1) 
with the number of spins n (cf. [l^ ). This means that the signal-to- noise in 
the spectra decreases exponentially with n, which precludes extending this em- 
ulation beyond 8—12 spins in the foreseeable future. Subject to these same 
limitations, however, an NMR computer is also able to directly estimate the 
expectation values of its observables . While this does not lead to any asymp- 
totic performance gains over what can be done with a quantum computer, it 
remains a potentially significant advantage. We conclude that the computa- 
tional potential of NMR spectroscopy, and of ensemble quantum computing 
more generally, has yet to be fully explored. 
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